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Q I Abstract: we study fermion masses and flavor mixing in a supersymmetric 5*0(10) 

^ ' model, where 10, 120 and 126 Higgs multiplets have Yukawa couplings with matter 

multiplets and give masses to quarks and leptons through the breaking chain of a 
Pati-Salam group. This brings about that, at the GUT energy scale, the lepton 
, mass matrices are related to the quark ones via several breaking parameters, and 

qj I the small neutrino masses arise from a Type II see-saw mechanism. When evolving 

i-^ ' renormalization group equations for the fermion mass matrices from the GUT scale 

to the electroweak scale, in a specific parameter scenario, we show that the model 

S^ . can elegantly accommodate all observed values of masses and mixing for the quarks 

and leptons, especially, it's predictions for the bi-large mixing in the leptonic sector 
are very well in agreement with the current neutrino experimental data. 
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I. Introduction 

The flavor problem about fermion masses and mixing has been a open question 
in particle physics [1 . In the Standard Model (SM), the fermion masses and mixing 
angles are completely arbitrary and neutrinos are massless, but recent experiments 
uniquely specify that neutrinos have small masses and bi-large leptonic mixing angles 
[2]. The present global analysis give at 90% C.L. [S] 

/\ml^ ^ 7.1 X 10"^ eV^ \/\ml^\ ^ (1.3 - 3.0) x 10"^ eV^ 

tan^ 9i2 ^ 0.4 , sin^ 2^23 > 0.92 , sin ^13 < 0.23 , (1) 

which are distinctly different from large masses and small mixing angels in quark 
sector. Incorporating the diverse values of fermion masses and mixing into a theory 
with the small number of parameters is a large challenge for theoretical particle 
physicists. There have been many ideas proposed for the purpose Q, among others, 
supersymmetric (SUSY) models based on the grand unified theory (GUT) such as 
SUSY 5*0(10) models are theoretically well-motivated extension of the SM and have 
received growing attention 0. 

SUSY 5*0(10) models provide the most natural framework for generation of neu- 
trino mass and realization of see-saw mechanism ^6 , in addition, the fermion mass hi- 
erarchies are also understood very well by appealing to additional family symmetries 
[7j. However, details of symmetry breaking chains and Higgs field structures plays a 
crucial role in specifying the model and determining the fermion masses and mixing 
angles [HIHI- One can generally obtain the correct quantum number for the SM par- 
ticle content and implement the see-saw mechanism through the breaking chain of 
either the Pati-Salam group SU{4:)c®SU{2)l®SU{2)r (O422) or SU{h)®U{l) (G51). 
The former breaks the left-right symmetry SU{i)c ® SU{2)l ® SU{2)r ® U{1)b-l 
(O3221) at the B-L energy scale and gives the right-handed neutrino masses by means 
of the renormalizable couplings of the fermion fields in 16 spinor representations with 
the Higgs multiplets in 126 representation, in which the i?-parity is automatically 
conserved at all energy scales, whereas the latter makes use of the nonrenormaliz- 
able couplings of the matter spinor fields with the 16 Higgs multiplets to generate 
neutrino masses. The different symmetry breaking chains give rise to different mass 
relations between quark sector and lepton sector. However, a successful GUT model 
is required to naturally account for the difference of masses and mixing between 
the quarks and the leptons. The minimal SUSY 5*0(10) model jH], where only one 
10 and one 126 Higgs multiplets have symmetric Yukawa couplings with matter 
multiplets, seems difficult to consistently incorporate the realistic neutrino oscilla- 
tion parameters. It was recently realized that there are possibilities for reproducing 
observed fermion masses and mixing by considering new Higgs multiplets |10j . 

In this works, we extend the minimal SUSY 5*0(10) model by introducing anti- 
symmetric couplings of one 120 Higgs multiplets with the matter fields. Supposing 
the G422 breaking chain, the left-right symmetry G3221 is broken by the O422 com- 
ponent (10, 1, 3) contained in the 126 Higgs multiplets and large Majorana masses 



are given to the right-handed neutrinos. Below the B-L energy scale, this model is 
described as the minimal supersymmetric see-saw standard model with superfield 
couphngs of the left-handed lepton doublets with the SU{2)l Higgs triplets, which 
are from the (10, 3, 1) component of 126 and can give small Majorana masses to 
the left-handed neutrinos. The two Higgs doublets in the minimal supersymmetric 
standard model (MSSM) are linear combinations of the SU{2)l doublet components 
from different SO{10) representations of Higgses, which all contribute to electroweak 
symmetry breaking and give Dirac masses to all the fermions. The effective Majo- 
rana masses for the light left-handed neutrinos are generated through the Type II 
see-saw mechanism ^J. Because of receiving contributions from different compo- 
nents of the 120 Higgs multiplets, the mass relation between the up-quark sector and 
the Dirac neutrino sector is different from the mass relation between the down-quark 
sector and the charged lepton sector. As we will see below, the characteristic struc- 
ture of Higgs fields and symmetry breaking superpotential can lead to fit completely 
all the observed data of masses and mixing for the quarks and the leptons. 

The remainder of this paper is organized as follows. In Section II we outline the 
model and derive the GUT relation among the fermion mass matrices, and then the 
renormalization group equations (RGEs) of the mass matrices are introduced. In 
Sec. HI, a detailed numerical analysis of the masses and mixing angles of the neutri- 
nos is given in a specific parameter scenario satisfying the experimental constraints. 
Sec. IV is devoted to conclusions. 



II. Model, Mass Matrices and Renormalization Group Evolution 

We consider the SUSY model based on the GUT gauge group 5*0(10). The 
matter fields \E'^^ in one 16 spinor representation of the 5*0(10) group contain all 
the quarks and leptons as well as the right-handed neutrino of each generation. 
Higgs fields which can couple to spinor fermions at the renormalizable level include 
only the H^^, H^"^^, and H^'^^ multiplets of the 10, 120, and 126 representations 
under 5*0(10), respectively. The gauge invariant Yukawa couphngs superpotential 
are such as 

Wsom = yl^^l'H^'^f + yJf^fH^^f + yl}'^l'H'''^f , (2) 

where i,j are generation indices. By virtue of the gauge symmetry, the Yukawa 
couplings, yjj' and yjf^, are symmetric 3x3 matrices, while yj^^ is an antisymmetric 
one. In our model, all the Yukawa couplings are assumed to be real to keep the 
number of free parameters minimum. 

We suppose that the GUT gauge symmetry descends to SU {3)c<^SU (2) l®U {1)y 
(O321) of the SM through two intermediate symmetries which are orderly the Pati- 
Salam subgroup G422 and the left-right symmetry subgroup 03221- At first, the 
50(10) symmetry is broken down to G3221 by developing non- vanishing vacuum ex- 
pectation values (VEVs) (54) + (45) of two Higgs multiplets in Higgs superpotential. 



The subsequent breaking to 6*321 is achieved by the (10, 1,3) component VEV of 
^126 jj^ eq.(2), which is a singlet under (7321. In the framework of G321, accordingly 
the Yukawa superpotential are rewritten as 



WG321 = Qi 
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where as usually the notations Q,L,u,e, etc. denote the chiral superfields for the 
quarks and leptons, and the subscripts of various Higgs superfields show that they 
respectively originated from corresponding representations of G422; the lowest in- 
dices u, d stand for belonging to the up-type and down-type SU{2)l Higgs doublet 
of G321, respectively. Note that a Clebsch-Gordon coefficient (—3) is generated in 
the lepton sectors when the components (15, 2, 2) of G422 are involved in the super- 
potential, which plays a crucial role in obtaining the Georgi-Jarlskog relations [T^ . 
The last second couplings will yield the left-handed neutrino Majorana masses after 
the neutral component of the SU{2)l Higgs triplet HW ^. develops a VEV as small 
as the neutrino masses measured in experiments. The last term gives the right- 
handed neutrino Majorana masses, moreover, the VEV (H}^q^^-,) should be close 

to the GUT energy scale so as to successfully implement the see-saw mechanism. At 
last, the electroweak symmetry breaking is accomplished by the electroweak scale 
VEVs of neutral components of these Higgs doublets in eq.(3), through which all 
the fermions acquire Dirac masses. 

From the point of view of phenomenology, omitting the heavy Higgs fields de- 
coupling from the low energy (here we assume that a doublet- doublet Higgs mass 
splitting is realized), the two light Higgs doublets in the MSSM, whose VEVs are 
Vu, Vd, are realistically identified with some appropriate linear combinations of the 
above up-type and down-type Higgs doublets, respectively. This corresponds that 
the up-type Higgs doublet VEVs are in proportion to Vu, while the down- type ones 
are relative to v^. To sum up, the various Higgs field VEVs in eq.(3) are described 
by 
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(4) 



where f„ = sin Pv,Vci = cos/3f, and v = VIA GeV is the VEV of the Higgs field 
in the SM, Mq is the GUT energy scale. The coefficients c'^/^iik = 1, ... ,4) arising 
from the mixing of the Higgs doublets are in generally complex, through which CP- 
violating phases are introduced into the fermion mass matrices. The parameters 
cliCr are positive, by which the mass scales of the left-handed and right-handed 
neutrino are related to the electroweak and GUT energy scale f , Mq- 

After the above symmetry breaking, the lagrangian relevant to the fermion 
masses is now obtained as follows 



--^mass = ULiMl'jURj + duMf-dRj + eLiM^jCRj + ULiM^-URj 
+ lfL^Ml^,^l^ + -ifR^M^jVR, + h.C. , 



(5) 



where the Dirac mass matrices for the quarks and leptons, as well as the Majorana 
mass matrices for the left-handed and right-handed neutrinos are, respectively, given 
by 
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(6) 

After integrating out the heavy right-handed neutrino fields decoupling from theory 
below the B-L breaking scale, the effective Majorana masses of the light left-handed 
neutrinos are realized by the type II see-saw mechanism as 
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(7) 



All the fermion mass matrices in eq.(6) are characterized by two energy scales f , M^, 
three positive breaking vacuum parameters tan/3, cl, cr, three real Yukawa coupling 
matrices yjj' , yl^^ , yjj^ , and eight complex mixing coefficients c'^,c'^, but these pa- 
rameters are not independent each other. 

Owing to only the three Yukawa coupling matrices, the fermion mass matrices 
in eq.(6) are, at the GUT scale, correlated each other. We define the following 
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The matrices M^ ^ denotes the symmetric and antisymmetric parts of the fermion 
Dirac mass matrices, respectively. The dimensionless parameters r;,r4,r5 are vir- 
tually some ratios of the Higgs doublet VEVs in eq.(4). By virtue of the special 
definitions, it is easy to derive the GUT relations among the quark and lepton mass 
matrices as follows 
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The symmetric parts of the fermion mass matrices are expressed as linear combina- 
tions of Mg and M|, while the antisymmetric parts of them are linearly dependent 
on only M^. Various coefficients of the linear combinations are determined by the 
seven parameters ri, . . . , rs, r/,, r/j, through which the lepton mass matrices are pre- 
dicted by the quark mass matrices. 

In the following discussion, we consider the model in a specific scenario. First, 
adopting a generation basis in which the up-type quark mass matrix M" is real 
and diagonal. Secondly, assuming that the down- type quark mass matrix M"^ is 
Hermitian. Namely, the quark mass matrices are described as 



M^ = diag K, ml m°) , M' = U'ckm diag {ml m°, m°) {U'ckm^ 



(10) 



where U(j^j^j is the CKM matrix with three mixing angles and one CP-phase in the 
quark sector ^2]; the superscript means that the quark masses and mixing are 
evaluated at the GUT scale (hereafter as such). As a result, it is very easy to know 
from eq.(6) that the relation cL = — 3cu , and the coefficients cj^,cj^ are real, 

while c 



(3) M 



are pure imaginary. Lastly, making a supposition that the coefficients 



Cu ,Cu are also real. To collect these together, we can draw a conclusion that 
all the charged fermion mass matrices in eq.(6), M'^,M'^,M'^, are Hermitian, while 
all the neutral ones, M", M^,M^, are real. The CP-violating sources come from 
only the antisymmetric parts in the down-type quark and charged lepton sectors. 
Furthermore, It can be inferred by eq.(8) that r4 = —3, and ri,r2,r^,rL,rFt are real 
(moreover, tl and rn have the same sign) but r^ is a pure imaginary. In terms of 
eq.(9), now this six parameters in addition to the ten values of the quark masses 
and mixing in eq.(lO) determine fully all the mass matrices in the lepton sector. 
On account of the Hermite matrix M"^ satisfying the constraint equations 

TrM" = m'^^ + m^ + m'l , DetM^ = ml-m\-m%, 

Tr {{M^f\ = {mlf + (mj)^ + {mlf , (11) 

the three parameters ri,r2,r5 can actually be solved out via inputting the charged 
lepton mass eigenvalues. Although one can find two sets of real solutions, whose 
parts involved in ri,r2 are the same except the parts r^ are different, the two M'^ 
matrices determined respectively by them are only a transpose each other. Without 
loss of generality, we only consider the solution whose r^ is relatively larger in the 
following analysis. Therefore, the three masses of the charged leptons fix completely 
the parameters ri,r2,r^, sequentially determine the charged lepton mass matrix by 
eq.(9). 

Now the left parameters rL,rji,r^ are restricted through fitting the masses and 
mixing in the neutrino sector. In terms of eq.(7), the effective left-handed neutrino 
Majorana mass matrix M^^ is a real symmetric, but only the two mass-squared dif- 
ferences are known from eq.(l), so we consider the following simultaneous equations 



ri _ ™0 , _0 , _0 T\^^.^\JL 



TrMgg = ml + m)^ + m''^ , DetM^^ = ml ■ m)^ ■ m 



3 ; 



Tr 



{Mt^y 



{m^f + {mlf + {mlf 



A'ml = {mlf - {m\f , AX2 = K)' - {<? ■ (12) 

Taking /\^m\^^ /\^ 711^2 i^z as input parameters, the values of m\,m2.,m\,ri,rR are 
solved out. The normal (or inverted) hierarchy corresponds to l\^m\2 > (or < 0). 
Note that four sets of real solutions can be found for eq.(12), but only two sets 
satisfy conditions that tl^tr have the same sign. In fact, the two sets of solutions 
are only contrary sign. Putting them into eq.(9) and using eq.(7), they actually give 
the same matrix M^^, so one can adopt either of them. 

To summarize the above scenario, there are in all sixteen independent input 
parameters evaluated at the GUT scale. They are six quark masses, three angles 
and one CP-phase in the quark CKM matrix, three charged lepton masses, two neu- 
trino mass-squared differences and one free parameter r^. Providing these parameter 
values, all the fermion mass matrices at the GUT scale are entirely determined. Fur- 
thermore, running the fermion mass matrices from the GUT scale to the electroweak 



scale according to the following introduced RGEs, we can compare our results with 
the experimental data of the fermion masses and mixing at the electroweak scale. 

Below the B-L scale, the right-handed neutrinos are decoupled by the see-saw 
mechanism. The model particle spectrum are identical to ones of the MSSM with 
the effective left-handed Majorana neutrinos. We introduce the Yukawa coupling 
squared matrices for the charged fermions as follows 



V V 

5, = i±i^^^(M^)tM^ 



(13) 



and then the one-loop close RGEs for these Yukawa coupling squared matrices, the 
effective neutrino mass matrix and the gauge coupling constants of G321 are given 

by mi 
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(17) 



where x = ^^{Q/Mg), h = (33/5, 1, —3), and / is a 3 x 3 unit matrix. Note that 
the eq.(14) and eq.(15) close on themselves, in addition, we neglect the running 
effects for the right-handed Majorana neutrinos since their mass scale is close to 
the GUT scale. For given tan/3, the previous calculated mass matrices and one 
unified gauge coupling constant are taken as the input values at the GUT energy 
scale Qgut = Mq, consequently, we can solve the above RGEs numerically and 
achieve the values of ^^(x)? Sf{x)y ^esix) ^^ the electroweak energy scale Qweak = 
Mz- The electroweak scale fermion mass eigenvalues are subsequently obtained by 
diagonalizing the Yukawa coupling squared matrices and the effective neutrino mass 



matrix as follows 

UuSuiXwWl = diag (mlixw), mKx^), m^^ixw)) , 

V 

UdSd{Xw)U^d = i diag {ml{xw),m1{xw),ml{xw)) , 

UeSe{Xw)Ul = ^ diag (m^(x^), mj(x»),m^(x^)) , 

UuM^ff{Xw)U^ = diag (mi(x«,), m2(x«;), "^3(Xt«)) , (18) 

where Xw = ^t^{Mz/Mg) denotes that the fermion masses and mixing are evaluated 
at the electroweak scale. Accordingly, the quark and lepton mixing matrices are 
given by [131 CS] 

UuU^d = U'ckaAx^) , UMt = U'cKMiXn.) diag (e^^S e^^^ O) , (19) 

where /3i, /32 are two Majorana phases in the lepton mixing matrix. Finally, the mix- 
ing angles and CP- violating phases in the unitary matrices ^c'Wa/(Xu.) are worked 
out by the standard parameterization in ref. [TH] . 



III. Numerical Results 

In this section, we present numerical results of our model. First of all, we fix the 
electroweak scale, the GUT scale and the unified gauge coupling constants such as 

Mz = 91.2 GeV, Mg = 3.1 x 10^*^ GeV, ai(0) = ^2(0) = 03(0) = 0.0408 . (20) 

According to eq.(14), three gauge coupling constants at the electroweak scale are 
found to be 

«i(x») -0.0168, «2(X«;)~ 0.0335, «3(x») ~ 0.1172 . (21) 

These are in accordance with the current experimental measures very well [TB] . 

Secondly, the Higgs VEV v is also fixed to f = 174 GeV, and the two rep- 
resentative values of tan/5 are taken as tan/3 = 10 and tan/3 = 30. In the case 
tan/? = 10(30) (Note, hereafter as such, the values in the parenthesis are corre- 
sponding to the case tan/3 = 30), we input the quark and charged lepton masses as 
well as the quark mixing angles and CP-violating phase at the GUT scale as follows 
(in GeV) 

m° =0.00106(0.00107), m° = 0.307 (0.312), m° = 136(143), 
m° = 0.00125 (0.00135), m° = 0.0265 (0.0286), m° = 1.02 (1.18), 
m° = 0.000324 (0.00035), mj = 0.0685 (0.074), m° = 1.17 (1.32), 

s?2 = 0.2229 (0.2229), s°3 = 0.0348 (0.0339), s% = 0.003 (0.003), 5^^ = 59° (59°) . 

(22) 



After using eq.(14) and eq.(15), the equivalent values of the above masses and mixing 
at the electroweak scale are solved out such that (in GeV) 



^5 



m„(x») ^ 0.00233, m,(xj ~ 0.677, mi(x^) ^ 181, 
rudixw) ^ 0.00439, m,(x^) ~ 0.0929, mb(x^) ~ 3.01, 
me (x«,) ~ 0.000487, m^(xj~ 0.103, m,(xj ~ 1.75, 
sUxn.) ~ 0.2229, 4(x«,) ^ 0.0412, s^x^) ^ 0.0036, df^ix..) ^ 59°. (23) 

They are completely consistent with the current status of the quark and charged 
lepton masses and the CKM matrix at the Mz scale |16j . 

Lastly, we take the input values of the neutrino mass-squared differences at the 
GUT scale and the parameter r^ as (remember that r^ is a pure imaginary) 

A^mj^ = 1.39 (1.50) x 10"^ eV^ A°m^2 = 3.9 (4.3) x 10"^ eV^ 

rg = 759i (1105i) . (24) 

Here, we only consider the case of the normal hierarchy, namely A^rrigg > 0. For 
the case of the inverted hierarchy, a similar analysis can also be performed. After 
the RGBs of eq.(14)-(16) running simultaneously, the output values of the masses 
and mixing for the effective left-handed neutrinos at the electroweak scale are found 
as follows (in eV) 

mi ^ 0.0084 (0.0050), ms ^ 0.0119 (0.0098), mg ^ 0.0463 (0.0459), 

Amj^ ^ 7.1 (7.1) X 10"^ Aml^ ^ 2.0 (2.0) x 10"^ 

sin^ 2^12 ^ 0.816 (0.816), sin^ 2^23 ~ 0.960 (0.948), sin^ 2^13 ^ 0.154 (0.082), 

tan2 0i2 ^ 0.40(0.40), sin^ig ^ 0.20(0.14), 
(5^3 ^ -0.04l7r (-0.0397r), pi ^ 1.0457r (0.0457r), p2 ~ 0.4957r (1.4957r) . (25) 

The above results are excellently in agreement with the recent neutrino oscillation 
experimental data in eq.(l). The three neutrino masses are all within experimental 
limits. The values of 1712 and m^, are close to 0.01 eV and 0.046 eV, respectively, 
but mi varies with tan/3 about range (0.005 — 0.008) eV. The two mass-squared 
differences and the three mixing angles are all at the center values of the allowed 
regions. The leptonic CP- violating phase 6[^ is approximate to — 0.047r (^ —7.2°), 
which is small very much in comparison with the quark one. One of the two Majo- 
rana phases is about vr (or 0), while the other is 7r/2 (or 37r/2) or so. In addition, 
the right-handed neutrino masses are obtained straightforward by diagonalizing the 
mass matrix M^ in eq.(9) such that (in GeV) 

Ml ^ 9.1 X 10^2 (1.1 X 10^^), M2 ^ 1.7 X 10^^ (3.4 x 10^^), 

M3 ^ 1.7 X 10^^ (4.8 X 10^^) . (26) 

They are close to the GUT scale except that the Mi value is slightly low for the 
case tan/? = 10, which are just expected in the previous sections. 
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FIG. 1. For tan/? = 10, the three leptonic mixing angles as functions of the 
parameter r^/i . The red and blue curves correspond to Am'^2 ~ 2.0 x 10~^eV^ and 
3.0 X 10~^eV^, respectively, and they all hold Am^2 ~ ^-^ ^ 10~^eV^. 

Now we illustrate the numerical analysis in detail. In Figs. 1 and 2 the three 
leptonic mixing angles sin^ 26ij {ij = 12, 23, 13) are shown as functions of the model 
parameter rs/i for the case tan/? = 10 and tan/3 = 30, respectively. In the case 
tan/3 = 10 (30), we fix the input A^ml^ = 1.39 (1.50) x 10"^ eV^ at the GUT scale 
so that Am2i ~ 7.1 x 10^^ eV^ is always held at the electroweak scale. If the value 
of A°m^2 at the GUT scale is taken as A'^mj^ = 3.9 (4.3) x lO'^ eV^, which leads 
2.0 X 10~^ eV^ at the electroweak scale, this case is depicted by the 

3 



to Am32 ^ 



red curves in each figure. If A^m'^2 = 5.9 (6.4) x 10 "^ eV 



the value of Am|2 is 



accordingly increased to Am32 ~ 3.0 x 10" eV , this case corresponds to the blue 
curves in each figure. The other input parameters are given by eq.(20) and (22). 

The graphs show explicitly the bi-large mixing characteristic of the lepton flavor 
mixing, namely, a very large solar angle sin^ 2^i2 ~ (0.7 — 0.9), a nearly maximal 
atmospheric angle sin^ 26'23 ~ (0.9 — 1), and a very small CHOOZ angle sin^ 2^13 ~ 
(0.05 — 0.2). As the value oi r^/i is increased in the displayed range, the atmospheric 
angle sin^ 2^23 is almost horizontal and always hold the near maximum. The solar 
and CHOOZ angles go up slowly, but they still lie in the experimental bounds. 
The best fit for the experimental data is at r^/i ~ 760 (1100) for the case tan/3 = 
10(30). In addition, when the electroweak scale Aml^ varies from 2.0 x 10~^eV^ 
to 3.0 X 10"'^ eV^, the plots only shift finely without changing the overall trends. 
Finally, it is worth to point out that r^/i varying has almost no effect on the values 
of Arri^2 and Am'^2 at the electroweak scale, which are principally determined by the 
GUT scale A^mfg and A^m'^2^ respectively. In short, our model predict naturally 
the bi-large mixing structure, moreover, the three mixing angles and two mass- 
squared differences of the neutrino oscillations are realized simultaneously for certain 
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FIG. 2. For tan/? = 30, the three leptonic mixing angles as functions of the 
parameter r^/i . The red and blue curves correspond to Am'^2 ~ 2.0 x 10~^eV^ and 
3.0 X 10~^eV^, respectively, and they all hold Aml^ ~ 7.1 x 10"^ eV^. 

parameter regions. 



IV. Conclusions 

In summary, we have discussed the fermion masses and flavor mixing within the 
framework of the SUSY S'O(IO) model. In the GUT model, we introduce the anti- 
symmetric couplings of the 120 Higgs with the matter fields besides the symmetric 
couplings of the 10 and 126 Higgs in the minimal SUSY S'O(IO) model. The new 
Yukawa couplings and the left-right symmetry breaking chain not only bring about 
that the quark and lepton mass matrices are closely related to each other, but also 
lead to elegant explanation for the difference between the large flavor mixing for the 
leptons and the small flavor mixing for the quarks. After the renormalization group 
evolution from the GUT scale to the electroweak scale, all the current experimental 
data for the fermion masses and mixing as well as the gauge couplings are repro- 
duced correctly in our model parameter scenario. In particular, the bi-large mixing 
structure for the leptons is naturally arisen in certain regions of the parameter space, 
moreover, the solar and atmospheric mass-squared differences for the neutrinos are 
also accommodated simultaneously. All the nontrivial results indicate clearly that 
the SUSY S'O(IO) GUT model with the 120 Higgs couplings is worth investigating 
deeply for understanding the mystery of the fermion flavor. 
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